Fractional-filling Mott domains in two dimensional optical superlattices 
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Ultracold bosons in optical superlattices are expected to exhibit fractional-filling insulating phases 
for sufficiently large repulsive interactions. On strictly ID systems, the exact mapping between hard- 
core bosons and free spinless fermions shows that any periodic modulation in the lattice parameters 
causes the presence of fractional-filling insulator domains. Here, we focus on two recently pro- 
posed realistic 2D structures where such mapping does not hold, i.e. the two-leg ladder and the 
trimerized kagome lattice. Based on a cell strong-coupling perturbation technique, we provide quan- 
titatively satisfactory phase diagrams for these structures, and give estimates for the occurrence of 
the fractional-filling insulator domains in terms of the inter-cell/intra-cell hopping amplitude ratio. 

PACS numbers: 05.30.Jp, 73.43.Nq, 03.75.Lm 74.81. Fa, 



Originally introduced as a toy model of liquid helium 
trapped in porous media pj, the Bose Hubbard model 
is nowadays routinely brought to experimental reality 
in terms of ultracold bosonic atoms trapped in optical 
or magnetic lattices The power of this physics tool 
as a virtually ideal realization of the theoretical model 
stands in the broad range of parameters and configura- 
tions attainable. For instance, the hopping amplitude 
of the bosons across the lattice sites can be controlled 
quite directly by varying the strength of the laser beams 
providing the optical confinement. Additional flexibility 
comes from the possibility of tuning — and even revers- 
ing — the interparticle interactions via Feshbach reso- 
nances 0- Further aspects of ideality are the fact that, to 
all practical purposes, the system is at zero temperature 
and isolated from the environment. Thus, for instance, 
the phononic excitations typical of realistic condensed 
matter lattice systems are ruled out. This impressive 
degree of control over the system parameters played a 
key role in the breakthrough experiment in which the 
superfluid-insulator quantumphase transition predicted 
for the Bose- Hubbard model [lj was actually observed for 
a gas of ultracold atoms trapped in an optical lattice £j . 

A great variety of optical lattices can be attained 
through a suitable choice of the number and setup of 
the laser beamsproviding the optic confinement [Bj]. So 
far, ID 6], 2D Q, and 3D Q Euclidan lattices have been 
realized, as well as ID superlattices |8| and quasiperiodic 
lattices @. The progress in optical trapping techniques 
prompted a great number of proposals for non trivial 
optical lattices, including periodic geometries [Icj . qua- 
sicrystals 0, 0] and 2D superlattices . 

The complex periodic geometry of optical superlat- 
tices, entailing a multi-band single-particle spectrum, al- 
lows for fractional filling insulating domains in the zero 
tem per ature phase diagram of superlattice BH models 
[l3L Il4l ITU ITq . It has been furthermore observed that 
such domains may exhibit an unusual loophole shape, as 
opposed to the roughly triangular shape of the customary 
Mott lobes 0|. So far, the theoretical investigations on 



the quantum phase transitions in superlattice BH models 
focused on ID systems, where the fractional filling insu- 
lating domains can be directly related to the inter-band 
gaps in the single particle spectrum ^tJj owing to the 
mapping between the hard-core boson regime and free 
spinless fermions 0,^1- O n non strictly ID structures 
such mapping does not apply, and therefore the phase 
diagram cannot be inferred from the single particle spec- 
trum. Qualitative results have been obtained in the ex- 
treme regime where the hopping amplitudes among dif- 
ferent cells are much smaller than those within the same 
cell 0| . Quantitative results can be obtained resort- 
ing to numerical methods such as quantum Monte Carlo 
|2Ct I2H or density matrix renormalization group algo- 
rithms [22]. These approaches necessarily address finite 
systems, whose size has to be sufficiently large to pro- 
vide a good approximation of the thermodynamic limit. 
On superlattices the dimensional scaling of the computa- 
tional demand is made even more serious by the fact that 
the building blocks are cells comprising several sites, as 
opposed to single sites. On the other hand, satisfactory 
quantitative results have been obtained for (/-dimensional 
regular lattices based on analytical |2^| and numeric [24[ 
strong-coupling perturbative expansions (SCPE). A sim- 
ilar technique, adapted to the complex periodicity of 
superlattices [25|, has been recently introduced for the 
study of the ID structures whose experimental realiza- 
tion is described in Ref. @. 

Here we extend this technique, called cell SCPE, to 
more complex superlattices. We focus on the two realis- 
tic structures sketched in Fig. ^ namely the trimerized 
kagome lattice (TKL) proposed in Ref. [l3| and a two-leg 
ladder (TLL), which in some sense represents the mini- 
mal diversion from strictly ID systems. Such a structure 
can be realistically obtained by superimposing a direc- 
tional confinement on a suitably chosen portion of a 2D 
optical lattice, as it was done in Ref. Q in the case of a 
ID lattice, or combining simple lattices with different lat- 
tice constants, as proposed in Ref. j2(J. We observe that 
the strong interaction regime of the ladder BH Hamilto- 
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FIG. 1: Cell structure of the TLL (left) and TKL lattice 
(right). Dark and light thick solid lines represent intra-cell 
and inter-cell hopping amplitudes, respectively. A light gray 
shade signals the nearest neighbors of the central cell, which 
has a darker background. This shows that cell-lattice of the 
TLL (TKL) is a ID euclidean lattice (2D triangular lattice). 



the first referring to the cell it belongs to, and the second 
denoting its position within such cell. Thus, the hopping 
amplitudes across sites belonging to adjacent ^-site cells 
can be described introducing a set of I x £ matrices t cc . 
The row index of such matrices refers to the cell denoted 
by the first superscript, whereas the column index refers 
to the cell relevant to the second superscript. Symmetry 
considerations lead to conclude that there are z/2 such 
matrices, and that t c c ' ''''' '■' 



For instance, in the case 



of the TLL, where z = 2, one has only t^ c+ = t5ij. A 
further £x£, symmetric, matrix Tjh allows the description 
of the hopping amplitudes among sites belonging to the 
same unit cell. Thus, the BH Hamiltonian on a generic 
superlattice reads 



h = h + t v, h q = J2h c 



(1) 



nian maps onto a spin ladder model, a system that has 
received wide attention in the literature [53, Hj| ■ This 
once again shows that ultracold atoms in optical lattices 
provide a direct realization of condensed matter models 
over a wide range of model parameters |2fi| | . 

As expected, we find that the structures we consider 
exhibit insulating domains at critical fractional fillings of 
the form k/t, where k is a positive integer and £ is the 
number of sites in the unit cell (£ = 2 for the TLL and 
£ = 3 for the TKL). Since we assume that there are no 
energy offsets between the lattice sites, the non integer 
fillings correspond to loophole domains Our aim is 
to investigate the conditions for the occurrence of such 
domains. Before getting into details, let us briefly list 
our main results. Unlike the ID case, on the TKL the 
superlattice structure is not sufficient for the occurrence 
of fractional-filling domains. That is to say, the inter-cell 
hopping amplitude has to be sufficiently smaller than the 
intra-cell amplitude for these insulating domains to oc- 
cur, despite the fact that an arbitrarily small difference 
between these amplitudes is sufficient for opening a gap 
in the single-particle spectrum of the TKL. Conversely, 
on the TLL our results suggest that fractional filling in- 
sulating phases occur even if the single particle spectrum 
has no gap. These discrepancies are consistent with the 
fact that on non strictly ID structures the hard core limit 
of the Bose-Hubbard model is not equivalent to a free 
spinless fermion model [l8| . 

Let us now introduce the notation for the superlattice 
BH Hamiltonian describing the optically trapped ultra- 
cold boson systems we are interested into. The spatial 
arrangement of the unit cells inherent in a superlattice 
is described by the so-called adjacency matrix A of the 
cell-lattice, where A cc i is 1 if the cells labeled c and d are 
adjacent, and zero otherwise. As it is clear from Fig. ^ 
the cell lattice is a ID chain for the TLL, and a trian- 
gular lattice for the TKL. The homogeneity of the cell 
lattice entails that the coordination of the generic cell c 
is actually cell independent, z c — Y^' c A cc > = z. A generic 
site of the superlattice can be labeled with two indices, 



where, introducing the the boson operators at the j th site 
of the c th cell, a c j, a^ rj and n CJ = a^ja c _j, 



He = E 



t 

j,h=l 

refers to an isolated cell, while 



V 



t jj' a tj ac 



(2) 



(3) 



takes into account the interaction among adjacent cells. 
The parameters U > 0, \i and Vj appearing in Eq. @ 
are the boson-boson (repulsive) interaction, the chemi- 
cal potential and the energy offset of the j th site within 
the unit cell (2i|. We observe that the partitioning of a 
superlattice into a homogeneous lattice of identical unit 
cells, and hence the arrangement of terms in Eq. IjTjI. is 
not unique. However, the perturbative approach we are 
interested into suggests that a convenient choice is such 
that the inter-cell amplitudes are smaller than the intra- 
cell amplitudes. 

The zero temperature phase diagram of the superlat- 
tice BH model in Eq.Q) is standardly obtained in terms 
of the ground state energies relevant to different fillings 
. We recall indeed that Hamiltonian H commutes with 
the total number operator, . n c j, and hence it can be 
studied within a fixed-number sector of the Fock space 
without loss of generality. Denoting by En the ground 
state energy of H relevant to a population of N bosons, 
the boundaries of the insulator domain possibly corre- 
sponding to a given filling / are standardly obtained as 
fj,± = ±(EfM±i — Ejm), where M is the number of sites 
in the lattice. The insulator domain actually exists if the 
inequality fi + > /i_ holds strictly 

As we mention above, we provide the insulator domain 
boundaries fi± in terms of strong coupling expansions for 
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the ground state energy of Hamiltonian QJ, assuming 
that the inter-cell hopping term in Eq. © is the pertur- 
bative quantity. As it is easily understood, the eigenvec- 
tors and eigenvalues of the (identical) cell Hamiltonians 
are the key quantities in our perturbative expansions. 
Denoting by \k, N) the fc' h eigenstate of the cell Hamil- 
tonian relevant to a population of N bosons and Ef? the 
corresponding eigenvalue, the generic eigenstate of the 
unperturbed Hamiltonian Hq relevant to a population of 
Y] N c bosons and the corresponding energy can be writ- 
ten as 

|k > N))=®|fc ,JV c ) CJ £7(k,N)=X;< e , (4) 

c c 

where the subscript c labels the cells. In particular, set- 
ting k c = 1 and N c = L, Eq. @ describes the un- 
perturbed ground state relevant to the fractional filling 
/ = L/£, where we recall that £ denotes the number of 
sites in each of the identical cells. As we mention above, 
in order to describe the insulator domain boundaries for 
such filling, we also need the ground state energies for 
the so called defect states j2j| , obtained increasing or de- 
creasing the total population by a single boson. Equation 
Q gives the correct unperturbed result for such states, 
provided that the population of one of the cells is in- 
creased or decreased by one boson. However, such energy 
is clearly degenerate, since the population variation can 
affect any cell. This means that the defect states must 
be treated according to degenerate perturbation theory. 
We carried out the expansions up to the second order 
for both the TLL and the TKL, obtaining the ground 
state energies for the fractional fillings / — Ljl and the 
relevant defect states in terms of the unperturbed cell 
energies and cell matrix elelements (k, N\ah\k', N'), 
where ah is a generic annihilation operator appearing in 
Eq. © and N, N' = L - 2, L - 1, L, L + 1, L + 2. The 
form of the perturbative terms is quite similar to those 
described in the appendix of Ref. [25j in the case of ID 
superlattices, although more complex, essentially due to 
the more complex topology of the structures in Fig. ^ 
Note indeed that the perturbative term in Eq. © is 
such that a boson operator at a given site of a cell may 
be involved in more than one inter-cell hopping, unlike 
the ID case. This (not excessive) increase in complexity 
is greatly rewarded by a very low computational effort, 
mostly required by the complete diagonalization of the 
cell Hamiltonian J5J for a few values of the cell popula- 
tion. Since the cells of the structures in Fig. n com P r i se 
just two and three sites, we are able to provide satisfac- 
tory quantitative results at a negligible computational 
cost also for insulator domains relevant to fillings that 
would be prohibitive for brute force numerical methods 
such as quantum Monte Carlo or density matrix renor- 
malization group simulations. More in general, the por- 
tion of phase diagram of a superlattice accessible to our 
perturbative technique depends only on the number of 
sites within a unit cell, and not by the topology or dimen- 
sionality of the cell lattice. We checked the correctness 
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FIG. 2: Half-filling insulator domain for a TLL with intra- 
cell hopping amplitudes r and inter-cell hopping amplitudes 
t/2. The 2 nd order cell SCPE result (solid lines) is compared 
to quantum Monte Carlo simulations for a TLL comprising 
50 rungs (errorbars) . 



of our perturbative results against brute force numeric 
results (Lanczos or quantum Monte Carlo algorithm) on 
small periodic ladder and kagome lattices (comprising up 
to 9 cells). Our formulas hold for the most generic super- 
lattice described by Eqs. This means that the 
intra-cell hopping amplitudes, as well as the energy off- 
sets Vk, may be different from each other. Furthermore, 
there can be z/2 different inter-cell hopping amplitudes, 
where z is the coordination number of the cell lattice. 
Here we assume that Vk = 0, so that the fractional fill- 
ing insulator domains have a loophole shape |17|. Fur- 
thermore we assume that all of the intra-cell hopping 
amplitudes are equal to r, and all of the inter-cell hop- 
ping amplitudes are equal to r'. In these conditions the 
presence of loophole domains can be related to a simple 
parameter, i.e. the hopping ratio t'/t, that is basically 
the perturbative quantity in our expansions. This ratio 
can be also connected to the energy gap in the single 
particle spectra of the structures under concern, that on 
ID structures bears a strict relation to the presence of 
loopholes, due to the exact mapping between hard-core 
bosons and free spinless fermions (I7lll8| . 

In the case of the TLL, the gap in the single particle 
spectum disappears for r/r' < 2, whereas for the TKL 
it is present for any r ^ r'. Detailed investigations at 
half- filling for the ladder and at 1/3- filling for the TKL 
show that the thresholds for the occurrence of the rele- 
vant loophole insulator domains differ from those of the 
single particle gaps. Actually, Fig. [21 shows that, ac- 
cording to both cell SCPE and quantum Monte Carlo 
simulations, the half-filling loophole domain of the lad- 
der is still present for r/r' = 2, when the single particle 
gap is not present any more. The comparison between 
these results also shows that cell SCPE provide quanti- 
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FIG. 3: Phase diagrams for the structures in Fig. Top: 
TLL with t = 3r'. Bottom: TKL with r = 6r'. The fillings 
of the insulator domains (gray) are also given in the plot. 



tatively satisfactory results also when r' is a significant 
fraction of r, i.e. when the perturbative term is far from 
being infinitesimal. Furthermore, both cell SCPE and 
quantum Monte Carlo simulations for hard-core bosons 
suggest that the half-filling loophole domain of the ladder 
persists for t/t' < 2, and disappears when the hopping 
ratio becomes smaller than a finite quantity smaller than 
2. This partially agrees with the results reported in Ref. 
[28| for a spin ladder that maps onto the hard-core limit 
of our model for a suitable parameter choice. 

On the TKL the situation is in some sense reversed, 
since the loophole domains disappear when a single par- 
ticle gap is still present. More precisely, the gap is open 
for any r' < t, but it must be r' < r/c for the loophole 



insulator domains to occur, where c = Ci ~ 2.32 and 
c = C2 ~ 2.87 according to first and second order per- 
turbative results, respectively. We mention that these 
values of c are the same for the first three loopholes, rel- 
evant to filling 1/3, 2/3 and 4/3. Note that the fact 
that C2 > ci is coherent with our perturbative approach, 
whose espansion parameter is basically r'/r. 

We conclude by commenting Fig. [21 which shows the 
phase diagrams for the structures under investigation as 
provided by our second order cell strong coupling expan- 
sions. The value of the (fixed) hopping ratio t/t' is 3 
for the TLL (upper panel) and 6 for the TKL (lower 
panel). Note that, since the hopping amplitude depends 
exponentially on the strength of the optical lattice, the 
above ratios correspond to a relatively small modulation 
of the potential barriers between neighbouring sites. We 
emphasize that these ratios correspond to an arbitrarily 
small modulation of the potential profile of the superlat- 
tice. Indeed, we recall that r ~ e~ , r' ~ e~ Ih , where 
/ is the strength of the optical lattice, while h and h' 
are scaling factors related to the height of the potential 
barrier between adjacent sites belonging to the same cell 
or to different cells, respectively Q- It is hence clear that 
a sufficiently large / produces the desired ratio, provided 
that h < h' . This suggests that the fractional filling insu- 
lating phase can be always reached if the lattice strength 
is increased while keeping the potential profile fixed. In 
this respect we emphasize that Figs. [5] and |31 correspond 
to a slightly different approach, since we assume that t/U 
is varied while keeping the ratio t/t' fixed. This can be 
done by adjusting the setup and intensity of the laser 
beams givin g ri se to the optical superlattice, at least in 
principle dill IS- 

However it may prove more con- 
venient to keep the beam configuration fixed (i.e. to fix 
t and t'), while varying the interaction strength U, e.g. 
via Feshbach resonances 0. 
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